Interplay between couplings and common noise in phase synchronization: 
disagreement between global analysis and local stability characterization 
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We consider two coupled phase oscillators in the presence of proportional ("common") and inde- 
pendent white noises. The global synchronization properties of the system are analytically studied 
via the Fokker-Planck equation. When the "effective coupling" is big compared to the common 
noises, the former favors and the latter hinder synchronization. On the contrary, when the coupling 
is small compared to the proportional noises, we find that the latter induce synchronization, opti- 
mally when their intensities are big and in the n:m synchronization ratio. Furthermore, in such case 
a small value of the coupling is better for synchronization. Finally, we show that synchronization, 
which is a global property, must not be studied via local stability such as with Lyapunov exponents. 
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Synchronization of interacting oscillatory processes is 
ubiquitous in nature For any general system of cou- 
pled oscillators, phase dynamics approach can be usefully 
applied provided that the intensities of all interactions on 
the oscillators, such as couplings between oscillators, ef- 
fect of noise on the oscillators, etc, are small compared to 
the natural frequencies In such case, synchronization 
as an adjustment of rhythms of phase coupled oscillatory 
processes because of the interaction with each other is 
well understood. Another more intriguing kind of syn- 
chronization is the one induced by noise, which is ob- 
served in many natural and experimental systems, such 
as lasers 0, neurons or ecological systems The 
first theoretical approach in the study of noise-induced 
synchronization in phase oscillators was by analyzing the 
Lyapunov exponent Q , but Lyapunov exponent only de- 
termines the local stability. A well established method 
of analyzing the global behavior of the system is via 
the Fokker-Planck equation 0], and synchronization in 
uncoupled phase oscillators with noise has already been 
studied using this approach [H, [t| . 

Still, an analysis that takes into account the two routes 
to synchronization - interactions between the oscillators 
(couplings), and common noise -, and studies the in- 
terplay between each other is lacking. In fact, in most 
real systems, coupling among the oscillators and the ac- 
tion of noise are simultaneously present, e.g. [ToL 11|. 
In this letter we systematically investigate for the first 
time the effect of couplings, independent noises, and pro- 
portional noises (i.e., common noises but with different 
intensities) on phase synchronization. In particular, we 
show analytically and numerically how coupling and com- 
mon noise compete in achieving synchronization: two dif- 
ferent routes to synchronization are possible depending 
on whether the couplings or the common noise prevail. 
Moreover, the properties of the system when synchro- 
nized via couplings or by common noise clearly differ 
- for example, whether increasing the couplings or the 



common noise results in stronger or weaker synchroniza- 
tion -, and therefore our results will be useful for taking 
apart systems in nature synchronized via interaction or 
via common noise. 

The global behavior of the system will be studied via 
Fokker-Planck. After reducing and solving analytically 
the equation, we will analyze the role of non-common 
noise, proportional noise, and couplings in synchroniza- 
tion, and we will show the competence between couplings 
and common noise. Finally, we will present a case that 
shows that the Lyapunov exponent is not always a good 
indicator of synchronization. 

Let us consider a system of two coupled oscillators in 
the presence of proportional and independent noises. The 
equations for the amplitude variables read: 



x a (t) = F a (x a (t)) + e a V a (xi(t),x 2 {t)) 

(*)) y/l£Z(t) + H a (x a {t)) y/H£ri a (t), 



(1) 



for a=l,2. Here x a is the amplitude component of the 
a-th oscillator, F a is its individual dynamics, the V's 
are the coupling functions, and the e's are the coupling 
intensities. £(f) is the common noise, and rj a (t) are the 
non-common noises. G and H represent the coupling of 
the oscillators to the noises. £(i) and rj a (t) are assumed 
to be independent Gaussian white noises, with zero mean 
value and unit intensity (as the actual intensities were 
taken aside into \/D a an d \/E a ): (£(f)£(s)) = 5(t — s), 
( Va (t) r) (s)) = S aP 5(t - s), (£(f) r] a (s)) = 0. 

As long as the intensities of the couplings and of the 
noises remain small, we can make phase dynamics. Ap- 
plying the phase reduction method to equation ([T]) we ob- 
tain the following equations for the phase variables @, @] : 

4> a (t) =uj a + e a U a (0i(f),0 2 (f)) 
+ \^Z(t) g a (Mt)) + \/E^V a {t)h a (<p a (t)). (2) 

Here oj a is the natural frequency of the a-th oscillator. 
Note that, for the amplitude equations, the noises may 
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be either multiplicative, or just additive by making some 
of the functions G and/or H equal to 1; but the noises 
will become multiplicative for the phase dynamics. 

We are interested in how much synchronization of a 
given order n:m jlj equations © yield. We usually quan- 
tify the phase synchronization by evaluating how close 
to a fixed constant value the generalized phase difference 
0_ stays over time. We thus make a biyective change of 
variables from fa and fa to fau and 0_: 



1 d 



m<pi - n(p 2 , 



<)+ = mcpi + n<p 2 , 



6+ + 0_)/(2m), h 



f>+-fa-)/(2n). (3) 



The equation for 0_ is obtained by substituting (|2|) into 
the first equation of |3]). The functions v a , g a and h a 
in © must be 27r-periodic in all the phases, and there- 
fore they can be written as a Fourier series. Now we 
will make a simplification, by replacing such functions 
by the lowest-order relevant contributions: these simpli- 
fied equations will still yield all the important results. 
V\((f)i(t), fa(t)) will be replaced by aisin0_, because as 
long as fa- stays close to a given constant, other har- 
monics correspond to fast oscillations [lj. V2 will be re- 
placed by a2 sin<^>_. The lowest resonant terms for gi(fa) 
and 52(02) are respectively b\ sm.{mfa) and & 2 sm(n02) 
as we will show later [l2|. Finally, the h a have no res- 
onant terms as we will show later, so we just put the 
first harmonic Cxsin(fa) for h\ and C2sin(</>2) for hi- 
With only the relevant terms, equations ([2]) read fa — 
uj\ + £isin</>_ + A\ sin(rnfa) + B± rji(t) sin^>x, and 
<j) 2 = UJ2 + £2 sin fa- -I- A 2 £(t) sm(nfa) + B 2 772 (t) sin fa; 
where e a — a a e a , A a — b a a/ D a , and B a = c a \J E a . 
Using these equations and the first equation of ([3]), we 
get the equation <j>- = cj_ + e_ sin 0_ + [to A\ sin(m0i ) — 
nA 2 sin(nfa)} £(£) + mBi 771 (t) sin^i - nB 2 T72 (*) sin fa, 
where w_ = muj\ — nuj 2 is the detuning, and the "effective 
coupling" e_ is defined in the same way. 

The former expression is an Ito stochastic differential 
equation, whose associated Fokker-Planck equation will 
be analyzed, considering <j) + and fa- as the independent 
variables. Things simplify much if we realise that <fi+(t) 
is a fast varying variable: in the absence of couplings and 
noises, (f> a (t) = uj a t, so 4>+{t) = uj + t. In our case with 
couplings and noises, some corrections have to be made to 
this expression; but, as long as the coupling intensities 
E a , A a , B a are small compared to the frequency uj a , 
such corrections will be small compared to the quickly 
growing term lo + t, so we can take the approximation 
4>+{t) ~ ^+t, and we may therefore integrate out fau 
for one period in the Fokker-Planck equation; because of 
that, only the terms of the Fokker-Planck equation that 
do not contain any derivative with respect to fau will 
be relevant. The equation for the probability density 
W(fau,fa-,t) is given by Q: 



dW 
~dt 



<:) 



2-{[mii sin{mfa) — n A2 sin(n<5i2)] 2 W 



[mBi sin fa } 2 W + [nB 2 sin fa] 2 W} 



, (4) 



where the suspension points stand for the terms that in- 
volve at least a derivative with respect to <f>+. We will 
study the cases where the intensity of at least one noise 
is enough to make the system fall into the stationary 
solution, for which dW/dt = 0. In that state, we can 
take the approximation that the probability density W 
is almost independent of the fast variable </>+, because 
as we mentioned before fau s» lo + t, so t rs <p+ /lu+ and 
therefore = dW/dt w uj + dW jdfau , from which we 
get dW jdfav- sa 0. As a result, we can work with the 
probability density of fa-, P(<p-), which results from the 
integration of W within a period of fau. From the sec- 
ond line of ((3J, we see that the system is 4nm7r-periodic 
in the variables fau and fa-. P(fa-) is then defined 
as P(fa.) = (l/4nmir) f* nm7T W(fau,fa-)dfau » W(fa-), 
where the time dependence was taken out of W because 
of the stationary condition, and the (approximate) inde- 
pendence of W of fay was used in the last step. 

From ([I]), replacing W by P (</>_), writing fa and fa as 
a function of fau and fa- , and integrating out 4>+ between 
and 4nTO7r, we get: 



d 

dfa 



— (u>- + £- sin d>-)P ■ 
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[(to 2 (A 2 + B 2 ) 



dP 



dfa 



E- sm 



fa-)W) 



+ n 2 {A 2 2 + B\) - 2nmA 1 A 2 cos fa-) P] } = ^- = 0, (5) 

as the integrals of the terms that involve at least a deriva- 
tive with respect to 4> + are equal to zero because of pe- 
riodicity. The stationary condition was used in the last 
equality. The last term "— 2nm A\A2 cos</>_" in the left 
hand side of (O is the contribution to the integral of the 
term "— 2 m Ai sm(mfa) nA2 s'm(nfa)" coming from the 
first squared bracket in ([4]). That is the reason why we 
previously claimed that the lowest-order resonant terms 
of gi{fa) and g2{fa) in ^ - the couplings of the oscilla- 
tors to the common noise - were respectively b\ sin(m^i) 
and 62 sin{nfa): the contribution of this crossed term of 
the squared difference to the integral is zero except for 
the to and n harmonics respectively (and integer multi- 
ples [Hj])- So other harmonics just add their u A 2v terms 
to the "P 2 " terms arising from the non-common noises - 
see {A\ + B 2 ) and {A 2 + P|) in -, and as a result the 
non-resonant harmonics of the couplings to the common 
noise play esentially the same role as the harmonics of 
the couplings to the non-common noises. 

Equation (JSJ) reads that the expression inside the 
top-level bracket in the left hand side is a constant, 
called the probability current S. If we call fi(<f>-) = 
™4(w_ + e_sin^_) and / 2 (</>-) = to 2 (A\ + B\) + 
n 2 {A\ + B 2 ) — 2nm A\ A 2 cos fa-, equation §5§ yields 

p'ifa-) + [h{fa-) + U(fa-)] p(fa-)/f 2 (fa-) =s/f 2 (fa-). 

This equation is easily solved by using the integrating 
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factor M(0_) = exp[/(/ 1 + f 2 ')/f 2 ] = f 2 eMJA/M, 
so the general solution of the differential equation is 



P(0_) 



S I exp[V(x)]dx + N 

^ , (6) 



M<f>-) exp[K(0_)] 



where V(x) = f fi/ f 2 - The lower limits </>o and xo for 
the former integrals are arbitrary, but once such values 
are chosen, they have to be kept the same. For conve- 
nience, we will set both to be equal to — ir. The two 
constants N and S in (|6|) are obtained by imposing two 
conditions. The first one is that P ((/>-) has to be 2w- 
periodic (this comes from the fact that the equation to 
solve (JSJ is 27r-periodic in <f>~). The second condition is 
the normalization, i.e., the integral of P(</>_) within one 
period has to be equal to 1. 

fi(x) and f 2 (x) are 27r-periodic functions. A primitive 
V(x) of /1//2 can be obtained analytically, valid inside 
the interval — it < x < it. Also, the definite integral Iy 
of /1//2 within one period can be calculated as V(tt) — 
V(— 7r). Such integral Iy turns out to be proportional 
to the mismatch cj_. Putting all things together, the 
integral that defines V(x) can always be obtained for any 
value of x, by writing as x = (2k — 1)tt + tp, with k integer 

and < <p < 2tt: V{x) = hi h = f^T^ hi h + 

ig::$: +9 h/h = Kiv+j:: +v Mf2 = *iv+v{-*+ 

ip) — V(— 7r). It is clear, as Iy is defined as the integral 
within one period (or from the former expression), that 



V(x + 2air) = V(x) + aly, a € Z. 



(7) 



A relationship between the two constants N and S 
in j6| will be obtained by imposing the periodicity con- 
dition P(4>- + 2ir) = P(4>-). To write the expression 
for P{4>- + 2tt), we first write the integral in ((6| as 
J^; +27r eMV(x)]dx = I P + f*- +27r exp[V(x)]dx = I P + 

exp(P/) J_ exp[F(a;)] dx, where Ip = exp[V r (x)] dx, 
and was used in the last step. The first factor in the 
denominator of ([6|) is periodic, + 2ir) = f%(<f)_). In 

order to write the second factor exp[V(0_ + 2tv)], (O is 
used. Then, by imposing that P(<j>- + 2tt) and P (</>-) 
be equal, we obtain a relationship between S and N: 
S = [exp(Iy) — 1] N / Ip. Note that, for ly=0, or equiva- 
lently, the mismatch lu^ equal to 0, the probability cur- 
rent S vanishes. Writing S as a function of TV in ([6|): 



P(0_) 



N 



/ 2 (0-)e^-) 



1 



- 1 



Ip 



i V W dx 



(8) 

where the constant N is the one that makes P(4>~) to be 
normalized within one period, J P(</>_) d(j>- = 1. 

Some analytical results can be derived for the case of 
zero frequency mismatch. Such results hold, at least ap- 
proximately, as long as the mistach remains small. For 




FIG. 1: Effect of non-common noise on synchronization. The 
line is the probability density function from the theoretical re- 
sult ((8} - the normalization constant N has been numerically 
obtained -; the points come from simulation. The plots are 
for vi = 1.025, v 2 = 0.23, n = 9, m = 2, ei = 0.1, £2 = -0.02 
- the coupling term is e 1 ( 2 ) sin0_ in the equation for (j>i(2), 
so ei and E2 will have opposed signs in general -, Ai = 0.1 
and Ai = 0.02. B x and B 2 are 0.1 and 0.02 in (a), 0.3 and 
0.06 in (b), 1 and 0.2 in (c). From the simulation of plot 
(a), <f>- stays within a small interval, the system never falls 
into the stationary solution due to the weakness of the noises. 
Therefore, for (a) the solution of the Fokker-Planck equation 
must be understood as an average over noise and over initial 
conditions. In the same way, the probability density obtained 
in (a) from simulation is only valid for the initial conditions 
for which the simulation was performed. 



w_ = 0, the expression for the probability density func- 
tion becomes much simpler: Po (4>~) — Ni[m 2 (Al+Bf) + 
n 2 {A 2 2 + Bl)- 2nmA 1 A 2 cos(f>-]( 2e -/ nmAlA ^- 1 1 where 
N\ is a constant independent of 0_, to be adjusted so 
that the probability density be normalized to 1. This 
expression can be written in a more convenient way: 



Po 



-) = m 



m 2 B\ 



n z B,2 +(mAi- nA 2 ) 



\2 



-AnmAiAo sin 2 — 
1 2 2 



(2£_/iimAiA 2 )-l 



(9) 



Stronger synchronization will take place when the term 
depending on <j>- (the last one inside the bracket) out- 
stands out of the ^--independent terms (the other 
terms): the first three terms inside the bracket should 
be as close to zero as possible in order to obtain stronger 
synchronization. Two conclusions can be obtained from 
this. In the first place, the bigger the absolute value of 
B\ and P2, the weaker the synchronization: indepen- 
dent noise is always harmful for synchronization. This 
property is plotted in Figure Q] (in the figures, the fre- 
quency mismatch is not exactly zero, but small). Sec- 
ondly, the third term will be closer to zero the smaller in 
absolute value m A\ — nA 2 be. In other words: the val- 
ues of the couplings of the proportional noise to the first 
and second oscillators should ideally be in the same n:m 
synchronization ratio in order to enhace the strenght of 
synchronization. This property can be seen on Figure [21 
for (b), m A\—nA 2 is equal to 0, and the corresponding 
probability density is the most localized. Because of this, 
we will assume from now on that A\ and A 2 both have 
the same sign. 
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FIG. 2: Influence of the ratio between Ay and A 2 in synchro- 
nization. 1/1 = 1.025, v 2 = 0.23, n = 9, m = 2, ei = 0.001, 
e 2 = -0.0002, Si = 0.1, B 2 = 0.02 and A 2 = 0.1. Ay is 0.1 
in (a), 0.45 in (b), 0.8 in (c). 
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FIG. 3: Influence of the intensities of the common noise in 
coupling-induced synchronization. v\ = 1.025, v 2 = 0.23, 
n = 9, m = 2, si = 0.3, e 2 = -0.05, By = 0.1, B 2 = 0.02. 
Ai and A 2 are 0.225 and 0.05 in (a); 0.45 and 0.1 in (b); and 
0.675 and 0.15 in (c) - note that mAy — n A 2 = in the three 
plots -. Here we have for (a) the same as in Figure QJ a). 



After stating the former constraint for Ay and A 2 in or- 
der to improve synchronization, we now study how their 
intensities affect. Here we have two different scenarios: 
from the exponent in @ we can study the competence 
between coupling and common-noise in synchronization. 
When the "effective coupling" e_ is big compared to 
AyA 2 , more specifically e_ > nm A\ A2/2, the expo- 
nent is positive, so the probability density will have its 
maximum at 0_ = 7T. If the product Ay A2 increases, 
the bracket in ([9]) will depend more strongly on 0_ - 
therefore the contribution of the bracket will be towards 
improving synchronization - but, on the other hand, the 
exponent will be smaller - therefore the contribution of 
the exponent will be towards flattening the probability 
density, and therefore towards hindering synchronization 
-. As dU) is more sensitive to the second effect (the ex- 
ponent) than to the first effect, we conclude that, in this 
case of coupling-induced synchronization, bigger values of 
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FIG. 4: Synchronization increases with the intensities of the 
common noise when induced by the latter instead of by cou- 
plings, vi = 1.025, v 2 = 0.23, n = 9, m = 2, ei = 0.01, 
£2 = -0.002, By = 0.1, B 2 = 0.02. Ay and A 2 are 0.45 and 
0.1 in (a); 0.9 and 0.2 in (b); 1.8 and 0.4 in (c) - note that 
mAi — nA 2 = in the three plots -. 



the noise (either common or non-common) hinder syn- 
chronization. This property can be seen on Figure [3l 

When the intensities Ay and A2 of the common noise 
are big enough compared to the effective coupling, more 
specifically when e_ < nm Ay A 2 /2, we are in a different 
regime: common-noise-induced synchronization. As the 
exponent in ([9]) is negative, the probability density will 
have a maximum at = 0. Different from the previous 
case, a big value of Ay A 2 is good now in the two ways: 
it makes the (/>_-dependent term inside the bracket out- 
stand out of the other terms, and it makes the exponent 
be bigger in absolute value. This is shown on Figure [4] 
(and by comparing Figure [21(c) with a)). Therefore, 
we can have very narrow peaks in the probability den- 
sity and, as a result, strong synchronization, for small 
values of the non-commmon noises, and big values of the 
couplings of the common noise to the oscillators, whose 
"generalized difference" m Ay — n A 2 is small in absolute 
value (nevertheless, Ay or A 2 cannot be arbitrarily big 
because the approximation <fi + w u> + t would not hold 
any more; furthermore, from the physical standpoint, we 
must recall that the phase dynamics approach is valid 
when the values of the couplings are not big). In this 
regime of common-noise-induced synchronization, a big- 
ger value of the effective coupling e_ will result into an 
exponent smaller in absolute value: small values for the 
couplings between the oscillators are better for common- 
noise-induced synchronization. 

Finally, we show that local stability, usually stud- 
ied by means of the Lyapunov exponents, does not al- 
ways caracterize the synchronization properties of the 
system. In this letter we will just present a counter- 
example, and a more thorough analysis will be pub- 
lished elsewhere. Let us consider the most favorable 
case for common-noise-induced synchronization, where 
By = B 2 = m Ay — n A 2 = uj- = 0. As we have seen 
before, the probability density will have a maximum at 
4>- = 0. Equation |5]) becomes d#_ (— £_ sin <j)_ P) + 
d^JnmAy A 2 sin 2 (0_/2)P] = d t P. The correspond- 
ing Ito stochastic differential equation is: d<f>- = 

e_ sin dt + \^2nm Ay A 2 sin 2 (^_/2) dw, where w 
is a Wiener process By using Ito's lemma 

[lH, we get: (dlog|</>_|) = (e_ sin (j>- /</>- — 
nmAyA 2 sin 2 (0_/2)/0 2 ) dt. By taking the limit 0_ — > 
0, we get the Lyapunov exponent at 0_ = 0: 



A = e_ - 



nmAyA 2 



(10) 



Therefore, in the case nmAy A 2 /A < e_ < nmAy A 2 /2, 
the global analysis we have made via the probability 
density yields that the common noise prevails over the 
couplings in inducing synchronization, and that the sys- 
tem synchronizes with 0_ staying close to most of the 
time. Nevertheless, the Lyapunov exponent is positive at 
0_ = 0. This shows that the local analysis made via the 
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Lyapunov exponent sometimes fails when characterizing 
synchronization, which is a global property. 

In conclusion, we have studied phase synchronization 
between two oscillators in the general case: with cou- 
plings between each other, and the action of propor- 
tional and independent noises. Two different routes to 
synchronization are possible: via couplings and via com- 
mon noise. When one induces synchronization, the other 
hinders it, thereby showing the competence of the two 
ways of synchronization. In particular, we have stated 
a surprising property: when synchronization is achieved 
via common noise, smaller values of the couplings (ide- 
ally no couplings) favor it. Finally, we have shown that 
synchronization, that is a global property, must not be 
studied by means of a local stability analysis such as via 
Lyapunov exponents. This work was supported by EU 
Project BRACCIA. 
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